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A new method for studying the many-body response functions is elabo- 
rated and first applied to the ^He longitudinal response. An integral trans- 
form of the response function is calculated from the bound-state-type equa- 
tions for several versions of the N-N force. The equations arc solved with 
the help of the hyperspherical expansion. The final-state interaction is com- 
pletely taken into account. The results are compared with the integral trans- 
form of the experimental response function for 250 MeV/c< q < 500 MeV/c. 
The difference amounts to 20-25% with the experimental response substan- 
tially exceeding the theoretical one in the low-energy region. 
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Many calculations testify to the fact that the conventional form of the 
nuclear charge density is inapplicable to the description of the elastic form 
factors of three- and four-nucleon nuclei at q > 2.5 fm~^ values. However 
the elastic scattering occurs with a quite low probability at such q values 
and some non-typical nucleon configurations including those where all the 
nucleons arc close together may make a substantial contribution. In this 
connection it seems important to test a form of the nuclear 4-current in the 
inelastic processes and to study the (e, e') response functions. This requires 
a proper account of the nuclear final-state interaction. 

In this paper we present a microscopical analysis of the ^He longitudinal 
response function Ri with a full account of the final-state interaction. The 
conventional approach calculates the response functions directly from the 
definition 



and it requires obtaining the whole set of the complicated final-state contin- 
uum wave functions ■0/- (Here e ~ a; — q^/{2AM) is the nuclear excitation 
energy.) Therefore some model approximations always were used in such a 
calculation. We apply a new approach [1-3] which enables us to avoid calcu- 
lating %ljf and thereby to obtain accurate results once the underlying nuclear 
dynamics is specified. 

Define the reduced transition operator and the response function 



o = [Gj(g)]-^o, R{q,uj) = [Gj(g)]-2i?(g,^). 

Here = - a;^ and [4] G^{Q) = [1 + Q"" /{AM'')]-y^G^{Q) where 



is the proton Sachs form factor. We calculate the integral transform of the 
response 



instead of the response itself. We use the conventional single-nucleon expres- 
sion for the charge density O. Then to a very good approximation one can 
disregard the a;-dependence of the O operator and use the expression 
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where p„ = — Rem- It has been shown [1,5] that a) can be calcu- 
lated by first solving for the localized solution to the following inhomogenous 
equation 

{H-Eo + a)^ = 0% (4) 

where as in Eq. (1) ip^ is the ground-state wave function and Eq is the 
ground-state energy. In terms of \l/ we have [1,5] 

^{q,a)=<i>\d%>-a-'R,i (5) 

where Rq is the elastic contribution to the response. 

The solution to Eq. (^) is much easier to obtain than the functions ipf 
entering Eq. (1). Indeed, in contrast to the latter functions there is no need 
to impose the complicated large- distance boundary conditions in order to get 
the solution. Only the condition that the solution vanishes at large distances 
is needed in solving Eq. (^. Therefore methods that are used in solving 
bound-state problems can be utilized here. In particular for many-body 
systems Monte-Carlo Green functions technique can be applied. 

We have two possible ways to connect our theoretical calculations with 
experimental measurements. One way [5] is to compare a) with the same 
quantity obtained from the experimental -R(g, e) using Eq. (2). Another way 
[1] is to consider Eq. (2) as the integral transform and invert to obtain the- 
oretical -R(g, a;) and then compare the responses themselves. In the present 
work we use the first approach. It is worth mentioning that there exists 
a generalization [1] of this method to exclusive reactions including 2 ^ N 
reactions induced by strong interaction. 

In this first calculation we use effective central N—N forces [6-8] which are 
supposed to act in the s-wave. We supplement them with a reahstic singlet 
p-wave N — N force [9]. Its contribution is 4% at the lowest q values and it is 
negligible at the highest q values. The contribution from the triplet p-wave 
force is believed to be of the same size. Only the components of the proton- 
proton Coulomb interaction which are diagonal in the isospin T = 1/2,3/2 
quantum numbers are retained in the calculation. Even these components 
change the results at most by 3% at the lowest q values considered. 

Under these assumptions on the nuclear dynamics, Eq. (^) is split into 
independent sets of equations with a given orbital momentum L, and isospin 
T of the system. It is convenient to calculate the right-hand sides of these 
equations in the following way. Since ipo has L = then only the components 
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~ YLMLiPn) from the expansion of exp(iqp„) from Eq. (|]) contribute to the 
problem for a given L value. Let q be directed along the z axis. Then only 
the components with = give non-zero contributions and hence only the 
components of ^ with Ml = are different from zero. We have 

(H-Eo + cr)^LT = Olt^o, (6) 
Olt = \T X T\ E[^^ + ^^^]jL(gPn)lLo(pJ, (7) 

n=l ^ '^p \1) ^ 

*L,r(g, or) =< ^l.tVlt^q >, (8) 

oo 

$(g,a)=47r^(2L + l) ^LT^q.o). (9) 

L=0 T=l/2,3/2 

The functions '^lt have the same spin S" = 1/2 as ■i/'o- 

We solve Eqs. by an expansion in the hyperspherical harmonics. 
Denote K the hyperspherical momentum and [/] the type of symmetry of 
the spatial components of the basis functions. We have developed a com- 
puter code to construct complete sets of the basis functions with arbitrary 
K,L,S,T and [/] values using the Raynal-Revai transformation [10]. The 
coefficients of this transformation are evaluated using the recurrent formula 
of the K K + 2 type [11]. Although the net number of the basis functions 
with the same L, S, T and [/] values grows linearly with K there exists 
a possibility indicated in Ref. 12 to specify the basis states in such a way 
that only two of them contribute to the problem in our case if one disregards 
the Coulomb interaction. It is because the nuclear forces that we use here 
only act in two N-N orbital states. Only such states are retained in our 
calculation. 

We make a comment concerning calculation of $ at small a values. It is 
necessary to avoid large cancellations in the right-hand side of Eq. (|^). This 
is achieved if one uses in Eq. for L = 0, T = 1/2 the same K^ax value as 
that at calculating ip^. Then the pole terms cancel exactly. 

There exists a test which enables us to check the calculation as a whole. 
Namely, the leading term of $(g,(T) at high a values behaves as a^^. This 
term can be calculated independently from the sum rule, 

/■oo _ — _ _ _ 

lim^^^a^{q,a) = R{q,uj)de = Y,m < MO^O\il)o > -Rei- (10) 
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This allows one to check the right-hand side of Eq. . Besides the correct- 
ness of the calculation the test allows one to verify whether at high a values 
the results are stable against increasing Kmax and L^ax- 

The calculations are performed at g = 250, 400 and 500 MeV/c. The 
required accuracy of the calculated $ is determined by the accuracy of $ 
extracted from experimental data. The latter is predominantly determined 
by systematic errors of the data and it is typically [13] about several per 
cent. We use Lmax = 10 while Kmax = 20 for L = and Kmax = 10 for other 
L values. We have verified that these values are sufficient for our results to 
converge within the experimental accuracy. 

Let us first discuss the relative contributions from various L and T values 
to the right-hand side of Eq. (^. The relative contributions of lower L 
values increase as q decreases in accordance with Eq. (|^. For example, at 
cr — s> oo the L = contribution dominates at g = 250 MeV/c contributing 
about 60% to the net sum while at g = 500 MeV/c the L = 2 contribution 
dominates with the L = component contributing only about 12 %. Besides, 
the relative L = contribution proves to increase as a decreases. At a = 1 
MeV and q = 500 MeV/c the L = component contributes more than 
70%. At (T = 1 MeV and q = 250 MeV/c it contributes more than 98%. 
The reason probably is as follows. At small a values the values of R with 
low excitation energies enter the integral from Eq. (2) with the highest 
weights. At low excitation energies the hyperspherical centrifugal barrier 
hinders particles from being inside the reaction zone. Such a barrier is absent 
when K = 0, which value exists only for the L = component. The effect 
is more pronounced at lower q values since the spectrum is shifted to lower 
energies. 

The T = 3/2 contribution proves to be suppressed by one or two orders 
relative to the T = 1/2 contribution at L = 0. The reason is that K = 
is forbidden for T = 3/2. Indeed, the K = basis state is symmetrical 
under particle interchange. But the symmetrical spatial components of the 
final-state wave functions appear at T = 1/2 only. For higher L values the 
T = 3/2 contributions are typically several times smaller than the T = 1/2 
contributions as well. This may be due to the fact that the spatially sym- 
metrical final-state components provide more inter-particle attraction that 
increases the amplitudes of the final-state wave functions inside the reaction 
zone. 

We compare our results with a) obtained from ^He experimental 
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responses [13] using Eq. (2). We note that at g > 550 MeV/c the longitudinal 
data of Ref. 13 become less accurate and, besides, the sums over the data 
greatly exceed the theoretical sum rule values (see Eq. (|TUp). We therefore 
only consider the lower q values. In order to perform the integration in Eq. 
(2) with a sufficient accuracy and in particular to estimate the contribution 
from the unavailable high-e tails of the spectra we approximated the spectra 
by the following analytical expressions: a{uj — oJthreshY^"^ in the low uj region 
^thresh < 1^ < t^i, S^=o ^n"^" ^hc rcgiou of the peak ui < u < u}2 and 
X]^=o Cni^"*-"'*'"''' in the region beyond the peak. The parameters a,bn,Cn,a 
and LV2 are chosen from the least-square procedure. Additional requirements 
of continuity of the fitting spectra and their first derivatives at ui and 002 
points are imposed. It turns out that a good description at a quite wide 
ranges UJ2 < uj < oJmax of the spectra beyond the peaks are provided with 
a single term ~ cu"". (In case, say, exponentially decreasing tail-terms the 
description is worse.) The best a values range between about 4 and 5 for all 
the q values considered. Similar a values were found [14] in the ^He case. We 
extrapolate the fits obtained beyond the Umax values in order to take into 
account the contributions from the unavailable tails of the spectra. These 
contributions prove to be quite small in our case. They reach the maxima at 
high a values where they are between 1 and 2%. 

Four versions of central N — N forces acting in the s-wave are used in 
our calculation. These include the S2 and S3 potentials from Ref. 6, the 
MT(I+III) potential from Ref. 7 and the EH potential from Ref. 8. The 
S2 potential and the MT(I+III) potential reproduce the N — N low-energy 
properties and s-wave N—N phases up to high energies. The S3 potential fits 
the low-energy data and yields nearly correct values for the binding energies 
and rms radii of ^He and ^He. The EH potential fits the s-wave N — N phases 
up to high energies as well but it does not reproduce properly the low-energy 
A^- data (see Ref. 6). 

The results of our calculation for g=250 MeV/c and 500 MeV/c are shown 
in Fig. 1 and 2 respectively. For g=400 MeV/c the results are similar. We 
plot the quantity [^theoriQ, o') — ^exp{(l, c")]/^ea;p(?, cr) against a. At the lowest 
a values this quantity amounts to -25% for the S2 potential and -21% for 
the MT(I+III) potential in case of g = 250 MeV/c and -(22-23)% for both 
potentials in case of g = 500 MeV/c. For the more phenomenological S3 
potential it equals -17% at g = 250 MeV/c and -24% at g = 500 MeV/c for 
the same a values. For the EH potential the corresponding numbers are -8% 
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and -20%. (The drawback of the EH force in reproducing the N — N data 
is pointed out above.) At higher a values the differences between ^exp and 
^theor monotonocally decrease and they do not exceed 6.5% in their absolute 
values for all the potentals as a goes to infinity. (At q = 250 MeV/c they are 
about 2% in cases of S2 and MT potentials. The differences in 6% in case of 
q = 500 MeV/c are perhaps connected in part with some systematic errors of 
the data at high q values, cf. above.) The results obtained indicate that the 
theoretical response functions are lower than the experimental ones at the 
low-energy wings. The deviation from experiment obtained may be caused 
by some deficiencies of the adopted conventional description of the nuclear 
states or, in principle, by some deficiencies of the conventional nuclear charge 
density operator. Further investigations are planned in this connection. In 
any case the results presented here show that a microscopic study of the 
nuclear (e, e') response functions with the final-state interaction fully taken 
into account is accessible and fruitful. 

We arc very indebted to C. Marchand for providing us with the exper- 
imental data. Helpful discussions with G. Do Dang and Yu. E. Pokrovsky 
are gratefully acknowledged. 
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CAPTION TO FIGURES 



FIG. 1. The relative deviation of the calculated transformant ^theoriQ, 
of the longitudinal ^He response function (see Eq. (2)) from that obtained 
from the experimental data [13]. 1-S2 potential, 2-MT(I+III) potential, 3-S3 
potential, 4-EH potential (see the text). The q value is 250 MeV/c. 



FIG. 2. Same as Fig. 1 but for q = 500 MeV/c. 
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